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A TENT PITCHING SCHEME 
MOTIVATED BY FRIEDRICHS THEORY 


JAY GOPALAKRISHNAN, PETER MONK, AND PAULINA SEPULVEDA 

Abstract. Certain Friedrichs systems can be posed on Hilbert spaces normed with a 
graph norm. Functions in such spaces arising from advective problems are found to 
have traces with a weak continuity property at points where the inflow and outflow 
boundaries meet. Motivated by this continuity property, an explicit space-time finite 
element scheme of the tent pitching type, with spaces that conform to the continuity 
property, is designed. Numerical results for a model one-dimensional wave propagation 
problem are presented. 


1. Introduction 

A commonly used approach for constructing numerical methods to solve time-dependent 
problems is based on the method of lines, where a discretization of all space derivatives is 
followed by a discretization of time derivatives. The resulting methods are called implicit 
or explicit depending on whether one can advance in time with or without solving a 
spatially global problem. The study in this paper targets a different class of methods 
referred to as locally implicit space-time hnite element methods, which advance in time 
using calculations that are local within space-time regions of simulation. Examples of 
such methods are provided by “tent pitching” schemes, which mesh the space-time region 
using tent-shaped subdomains and advance in time by varying amounts at different points 
in space. 

Ideas to advance a numerical solution in time by local operations in space time regions 
were explored even as early as [20]. Recurrence relations on multiple slabs of rectan¬ 
gular space-time elements were considered in [T3|, whose ideas were generalized to non- 
rectangular space-time elements for beams and plates in |2|. These works are not so 
related to the current work as some of the more modern references. Closest in ancestry 
to the method we shall consider is found in [22] where it was called explicit space-time 
elements. The space-time discontinuous Galerkin (SDG) method was announced almost 
at the same time in im and continues to see active development jlHl |2Tl |28] . Against 
this backdrop, we highlight two papers that brought tent pitching ideas into the numeri¬ 
cal analysis community [101 [Ej- The questions we choose to ask in this work have been 
heavily influenced by these two works. We should note that the name “tent pitching” 
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Figure 1. Tent pitching (read column by column) 


has been traditionally used for meshing schemes that advance a space-time front [HIES], 
but in this paper tent pitching refers to the discretization scheme together with all the 
required meshing. 

To give an overview of what is involved in a tent pitching scheme, consider the case of 
a hyperbolic problem posed in one space dimension with time as the second dimension. 
Given a spatial mesh, we pitch a tent by erecting a tent pole (vertically in time) at a vertex, 
in Figure (Precise dehnitions of “tents” etc. appear later - see Dehnition 4. 


as 


In the 

plots of Figure [T| the horizontal and vertical dimensions are space and time, respectively. 
The height of the tent pole must be chosen small enough in relation to the hyperbolic 
propagation speed, so that the domain of dependence of all points in the tent remains 
within the tent’s footprint. We then use the given initial data to solve, by some numerical 
scheme, the hyperbolic problem restricted to the tent. Proceeding to the next vertex 
where the second tent is pitched in Figure we hnd that the initial data combined with 
the solution in the previous tent, provides inflow data to solve the hyperbolic problem 
there. Solution on the newer tents proceeds similarly. This shows the sense in which tent 
pitching schemes are locally implicit: they only involve solving local problems tent by 
tent. 

Having explained tent pitching schemes in general, we should now emphasize that the 
main result of this paper is not a new tent pitching scheme (although one is included to 
show relevance). Rather, this paper is mainly concerned with answering a few theoretical 
questions motivated by tent pitching schemes. Indeed, our main result is a characterization 
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of traces of a Friedrichs space on a tent-shaped domain and builds on the recent advances 
in Friedrichs theory [niziiniiii]. To explain the Friedrichs connection, we should hrst note 
that all the previous tent pitching schemes use non-conforming space-time discontinuous 
Galerkin discretizations. Design of tent pitching methods within a conforming setting, 
while holding the promise of locally adaptive time marching with fewer unknowns, pose 
interesting questions: What is the weak formulation that the tent pitching scheme should 
conform to? What are the spaces? What are the hnite element subspaces one should 
use? These questions form the motivation for this study and while attempting to answer 
them, Friedrichs spaces and their traces appear naturally, as we shall see. While we are 
far from answering the above questions for a general Friedrichs system, our modest aim in 
this paper is to provide some answers for a few simple problems in one space dimension. 

Accordingly, there are two parts to this paper. The hrst and the main part of the 
paper consists of Sections |2]-[^ While results of Sections and are applicable to any 
abstract Friedrichs system. Section focusses mainly on an advection example and its 
implications for hyperbolic systems. This leads to observations on the traces of certain 
Friedrichs spaces. The theory clarihes a weak continuity property of the traces at the 
points where inhow and outhow part of boundaries (dehned precisely later) meet. It is 
relevant in the tent pitching context because in the tent-shaped domains used in tent 
pitching schemes, inhow and outhow boundaries always meet. The second part of the 
paper, consisting of Sections and designs an explicit space-time hnite element scheme 
of the tent pitching type using the spaces and weak formulations motivated by the hrst 
part. The method we construct is a low order method that works on unstructured grids. 
On uniform grids, comparison with a standard low order hnite diherence method does 
not reveal any striking advantages for the new method, as we will see in Section Yet 
we hope that this study will pave the way to a better understanding of conforming tent 
pitching discretizations, the spaces involved, and eventually lead to high order methods 
on unstructured grids for multidimensional problems. We begin with some preliminaries 
on Friedrichs systems in the next section. 


2. Friedrichs systems 


Our approach is inhuenced by the modern take on the classical work of Friedrichs m, 
as presented in misiini. Let L be a Hilbert space over M with inner product (■, ■)l and 
norm || ■ ||i, and let D be a dense subspace of L. Suppose A and A are linear maps from 
H into L satisfying 

(2.1a) (A(/),'0 )l = (0, AV’)^, 

(2.1b) 3c>0: ||(A + i)0|U <c||0|U, V0 G H. 


Let Wq be the completion of H in the norm ||0||w = (||0||i + ||A0||) 


p/2 


Then, with L 


as a pivot Hilbert space, identihed to be the same as its dual L', we have 2) C Wq Y L = 
L' C Wq. It is now standard to extend A and A as bounded linear operators from Wq into 
L, i.e.. A, A G JZ{Wq,L). Moreover, it is also well-known that Assumption (2.1) implies 
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that they can be further extended to A, A G ^(-h, Wq) via 

(2.2) {AE,wo)wo — R’o)wo = (^, Vf' E L, wq ^ Wq. 

Here and throughout, we use {■,-)x to denote the duality pairing in X. Next, dehning 
W = {v E L : Av E L}, we observe that Wq C W and that W normed with || ■ ||vk dehned 
above is a Hilbert space. Hereon, the norm on any normed linear space X will be denoted 

by II ■ llv- 

An important ingredient in Friedrichs theory is the “boundary” operator D in L{W, W) 
dehned by 

(2.3) {Du,v)w = {Au,v)l — {u, Av)l \fu,vEW. 

This is an abstraction of an integration by parts identity. For any operator B E T(H^, W), 
we dehne B* E £(hF, W) by {B*u, v)w = {Bv, u)w for all u^v eW . For subspaces S <EW 
and R C W\ dehne the right and left annihilators by 

= {w' E W : {w', s)w = 0 for all s E S'}, 

= {w eW (s', w)w = 0 for all s' E R}. 

These results are well known [U]: 


Proposition 2.1. 


The following are consequences of Assumption (2.1).- 


(a) D* = D. 

(h) The norm ||ry||^p = (||w||| + ||y4tc|||)^''^ is equivalent to ||tc||w for all w E W. 

(c) kerD = IFq- 

(d) ranT) = W^. 


We will henceforth tacitly assume (2.1) throughout this section. In the traditional 
Friedrichs theory, another “boundary operator” M, also in L{W, W), plays a leading role. 
This is a generalization of certain matrices used by Friedrichs m to impose boundary 
conditions. In the generalization of Friedrichs theory to the Hilbert space setting, as 
described in [9], the operator M is assumed to satisfy 


(2.4a) {Mw,w)w^0, Vw G IT, 

(2.4b) IF = ker(T)-M)+ker(T) + M). 


The theory in |9] addresses the unique solvability of two problems: The hrst is to hnd 
a M G IF, given any / G L, satisfying Au = f (typically a partial diherential equation), 
and {D — M)u = 0 (typically a boundary condition). The second problem is the “dual” 
problem of solving Au = f satisfying [D + M*)u = 0. These two problems are uniquely 
solvable if and only if the following two conditions hold, respectively: 


(2.5a) A : ker{D — M) —>■ L is a bijection, and 

(2.5b) A : ker{D + M*) -E L is a. bijection. 

Some sufficient conditions for ( |2.5 ) to hold can be found in [Hi- 
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In pj, an intrinsic approach without the operator M was discovered. It uses the double 
cones 

= {w G hh : {Dw,w)w > 0}, 

C~ = {w E W : {Dw,w)w < 0}. 


The intrinsic approach replaces (2.4) by the following assumption on two subspaces of W 
denoted by V and V*: 

(2.6a) 1/ C C'+, 1/* C C", 

(2.6b) V=^D{V*), V*=^D{V). 


Clearly, (2.6b) implies that both V and V* are closed, and moreover, 
(2.7) keiD = WoCVnV*. 


Note that the reflexivity of Hilbert spaces and (2.6b) imply that 
(2.8) (VY = = D{V), 

1/^ = (^D(l/*))^ = D{V*). 

The theory in j9] provides sufficient conditions for unique solvability of two problems: 
The first is to find a u G IT, given any f E L, satisfying 

(2.9a) Au = f (typically a partial differential equation), 

(2.9b) u E V (typically a boundary condition). 

The second is the “dual” problem of solving for a. u E V* satisfying Au = f. These two 
problems are uniquely solvable if and only if 

(2.10a) A : V ^ L is a bijection, and 

(2.10b) A ■. V* ^ L is a bijection. 


A coercivity condition on A + A is sufficient for (2.10) to hold, as proved in [9]. However, 
for operators like the transient wave operator considered later, A + A is zero and cannot 
be coercive. 

Hence, our first point of departure from P] is the introduction of another simple suf¬ 
ficient condition for unique solvability. It requires that the operator A be bounded from 
below, a condition which is often easy to verify for time-dependent problems (see e.g. pU]). 
Although the next theorem requires both A and A to be bounded below, one of these 


conditions can be easily removed in most applications, as detailed in Remark 2.3 


Theorem 2.2. Suppose (2.1) and (2.6) hold. If there is a constant c > 0 such that 
A : V ^ L satisfies 

(2.11a) ^ c||m||w Vm G V, and 

(2.11b) ^ c||m||w Vn G V*, 


then (2.10) holds. 










6 


GOPALAKRISHNAN, MONK, AND SEPULVEDA 


Proof. Inequality (2.11a) implies that A-.V^Lis injective and has closed range. Hence 
H is a bijection if its adjoint is injective, i.e., if 


(2.12) [leL- (Hn,= 0 for all v eV} = {0}. 


To prove (2.12), consider an £ satisfying {Av.,tjL = 0 for all v ^ V. Then, for all 
wq G Wo C V, we have by (2.2) that {Ai, wo)wo = 0, from which it follows, by the density 
of Wo in L, that A£ = 0 and £ ^W. Hence we may apply (2.3), which yields 


(2.13) 


{Dv, £)w = {Av, £)l - (n, A£)l = 0, Vn e V. 


Thus £ G P(V) = V*, and so (2.11b) implies £ = 0. This proves (2.12). 
That H is a bijection is proved similarly. 


□ 


Remark 2.3. Note that under the assumptions of Theorem 2.2, if (2.11a) holds and A is 


injective, then (2.11b) holds with the same constant c. This is most easily seen by viewing 


H as a closed (possibly unbounded) operator on L with dom(H) = V. From (2.7), we 
know that T> C Wo G V C L. Since D is dense in L, the domain of A is dense in L. 
Hence the adjoint A! is a well dehned closed operator on L satisfying (Hn, s)/, = (n, H's) 
for all V G dom(H) and s G dom(y4'). By dehnition, dom(H') consists of all s G L for 
which there exists an £ G L with the property (s, Av)l = {£, v)l for all v G dom(y4) = V. 
In particular, whenever s G dom(H'), by (2.2), {s,Awo)l = {As,Wo)wo — (^)'?^o)l for 
all Wq G Wo G V, so As = £ and consequently s G W. Thus, whenever s G dom(y4'), 
both {s.,Av)l and {As,v)l coincide with {£,v)l for all n G H, which by ( 2.3| ), implies 
that {Dv,s)w = 0, which in turn implies that s G = ^D{V) because of (2.6|), i.e.. 


dom(y4') C V*. Combining with the easily provable reverse inclusion, dom(H') = V*. 
Next, we claim that A' is the same as A. Indeed, (Hn, s)l — (u, A s)l = (Tlw, s)vi/ = 0 for 
all V G dom(H) = V and s G dom(H') = V* = '^D(V), thus proving the claim. Now, 


since (2.11a) implies that the ran(H) is closed, by the Closed Range Theorem for closed 


operators P!. we conclude that ran(H') = ran(y4) is closed. Hence if A is also injective. 


then by standard arguments, (2.11b) follows 


To summarize, we have discussed two known approaches to abstract Friedrichs systems 
and introduced a new sufficient condition for unique solvability of Friedrichs problems. 


The hrst approach via (2.4) is closer to the classical theory (the M-approach) while the 


second is the approach via (2.6) (the H-approach). Whether these two approaches are 
equivalent is a natural question. It was shown in |9] that if an operator M exists that 


satishes (2.4), then V = ker(Zl — M) and V* = ker(Zl + M*) satishes (2.6). The converse 
remained unknown until it was proven in [I]. In the remainder of this paper, we will use 
only the H-approach. 
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3. A WEAK FORMULATION WITH BOUNDARY FLUXES 


Consider the following abstract boundary value problem: Given / G L and g G W, find 
u E W satisfying 


(3.1a) 

(3.1b) 


Au = f, 
u- g eV. 


Space-time Friedrichs systems with non-homogeneous conditions on space-time bound¬ 
aries (which includes initial conditions) can be abstracted into this form. 

To derive a weak formulation, we multiply (3.1a) by a test function v E W and use (2.3), 
to obtain {u,Av)l -|- {Du,v)w = This implies 

(3.2) {u,Av)l +{D{u-g),v)w = F{v), 
where 

(3.3) F{v) = (/, v)l - {Dg, v)w 


Now, we let D{u — g) in (3.2) be an independent “flux” variable q. This leads us to 
formulate the following variational problem: 


(3.4) 


Find u E L and q E (C*)"'' such that 
{u,Av)l + {q,v)w = F{v), 


Vn G W. 


The bilinear form on the left hand side will be denoted by b{{u,q),v). Our approach to 
the construction and analysis of this weak formulation is close (but not identical) to the 
approach in [3]. 

A similar derivation for the adjoint problem of hnding a. u E W, given f E L and 
g E VF, such that 


(3.5a) 

(3.5b) 


Au = /, 
u-gEV\ 


(3.6) 


Vn G W, 


suggests the following dual weak formulation: 

Find u E L and q E V~^ such that 

{u,Av)l - {q,v)w = F{v), 

where 

(3.7) F{v) = (/, v)l + {Dg, v)w 

The bilinear form on the left hand side will now be denoted by b{{u, q),v). 

In applications, the {■,-)w terms can typically be identihed as boundary terms, so q 
and q can be interpreted as boundary fluxes. Finally, note that by virtue of (2.8), we can 
equivalently use D{V) and D{V*) as the flux spaces in (3.4) and (3.6), respectively. 
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3.1. Wellposedness. Next, we prove that the new weak formulation is well posed and 
is equivalent to the classical formulation (3.1) in the following sense. 


Theorem 3.1. Suppose (2.6) and (2.10) hold. Then the following statements hold: 


(a) Given any F G W, there is a unique {u,q) G L x (1/*) that solves (3.4). Moreover, 
if F is as in (3.3) for some given f ^ L and g G W, then the solution {u,q) of (3.4) 
satisfies 


(3.8) 


Au = f, u-geV, q = D{u-g). 


(b) Given any F G W, there is a unique {u,q) G L x V-^ that solves (3.6). Moreover, if 
F is as in (3.7) for some given f E L and g G W, then the solution {u,q) of (3.6) 
satisfies 

Au = f, u-geV*, q = D{u-g). 

To prove this theorem, we will verify a uniqueness and an inf-sup condition in the 
following lemmas. 


Lemma 3.2 (Uniqueness). Suppose (2.6) and (2.10a) hold. Then, whenever u G L and 
q G satisfies h{{u, q),v) = 0 for all v G W, we have {u, q) = 0. 


Proof. Suppose 
(3.9) 


{u, Av)l-\-{< i,v)w = ^ Vn G lU. 


Since q G (U*) , we have (g, v)w = 0 for all v G Wq due to (2.7). Hence, choosing v = Vq E 
Wq in (3.9), we conclude that {u,Avo)l = 0. Hence, using (2.2), we have {Au,Vq)wo — 0 
for all Vq G Wq, which implies, by density, that Au = 0 in L. In particular, this shows that 
u is in W. We may therefore apply (2.3) to (3.9) to get {Au, v)l — {Du, v)w + {q, 'v)w = 0, 
for all V eW. Since Au = 0, 

{Dv,u)w = {q,v)w Vn G W. 

Since q G (U*)"'', the right hand side vanishes for all n G U*, so m G '^D{V*). Hence 
by assumption (2.6), u E V. By (2.10a), u = 0. Using this in (3.9), it also follows that 
g = 0. □ 


Lemma 3.3 (Inf-sup condition). Suppose (2.6) and (2.10a) hold. Then, there is a C > 
such that for all v E W, 

\b{{u,q),v)\ 


C\\v\\w < sup 

{u,q)eLx{V*A II 


Wv E W. 


Proof. By (2.10a), there is a c > 0 such that given any v E W, there is a unique w E W 
satisfying 


(3.10a) 

(3.10b) 

(3.10c) 


Aw = V, 
w eV, 

lu’Ww < cllnlU. 
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Then, since w & V = we have {Dv*,w)w = 0 for any v* G V*. Therefore, 

q = Dw is in (f/*)"*". Moreover, 

\\v\\l + \\Av\\l = {Aw, v)l + {Av, Av)l by ( |3.10aD 

= {w + Av, Av)l + {Dw,v)w by (2.3) 

Note that the right hand side eqnals b{{w + Av,q),v) as q = Dw. Continning, 

b{{w + Av,q),v) 


\vrLA\\Av\\i = 


< 


{w + Av,q)\\LxW' 

\b{{z,r),v)\ 


snp 

Az,r)&Lx(V*y 


{Z,r)\\Ly.W' 


w + Av,q)\\LxW' 

[w + Av,q)\\LxW'- 


By Proposition 


2.1 


b|) and ( |3.10c[ ), ||w + Av\\l < IIw^IIl + ll'^lliy < C'||n||^y for a C > 0 

□ 


depending on c. Moreover, if d denotes the norm of D, then 
Using these estimates to bonnd ||(w + Av, q)\\LxWA the lemma is proved 

Proof of Theorem 3jJ_. Lemmas 3.2| and 3^ verify the conditions of the Babuska-Brezzi 
theory, from which the stated nnique solvability follows. 

Now snppose F is expressed in terms of / and as in (3.3). Then choosing v = vo E Wo 
within the weak formulation, 

{u, Avo)l + {q, Vo)w = if, Vo)l - {Dg, Vo)w = (/, Vo)l 

we obtain {u,Avo)l = {Au,Vo)w = {fWo)L- This proves, by density, that Au = f in L, 
and consequently u G W. Then, returning to (3.4) and using (2.3) together with Au = f, 
we obtain {q,v)w = {D{u — g),v)w for all v G W, i.e., q = D{u — g). Finally to show 
that u — g E V = ~^D{V*), consider an arbitrary v* E V*. Then note that q E (U*)"*", so 

0 = {q, v*)w = {D{u - g),v*)w = {Dv*, u - g)w 


This proves (3.8). The remaining statements are proved similarly using (2.10b) in place 
of (2.10a). □ 


4. Examples 


The assumptions on which the previous theory is based can be verified for several 
examples. We begin with the simplest example in one space dimension in § 4.1 where 
all the ideas are transparent. We then generalize to the example of multidimensional 
advection in § |4.2| and establish a new trace theorem for the associated graph space. 


The hnal example in § 4.3 considers a general symmetric hyperbolic system in one space 


dimension and leads into the discussion on the wave equation in the subsequent section. 


4.1. An example with no space derivatives. We begin with a simple example in one 
space dimension that illustrates the essential points. Let K denote the open triangle in 
space-time {x,t) G M x M, with vertices at {x,t) = (0,0), (1,0), and (1,1). Set 

(4.1) L = L\K), ‘D = T){K), Au= — 
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(where Ti{K) denotes the set of compactly supported inhnitely differentiable functions 
on K). Obviously, A = —dt, so (2.1) is satished. We split the boundary of K into an 
inflow, outflow, and a characteristic part: 

diK = {(a;,t) G dK : t = 0}, doK = {{x,t) ^ dK : x = t}, 
dcK = {{x,t) e dK : X = 1}. 

Because dist{diK,doK) = 0, although the operator D is dehned on all W, we must 
be careful in speaking of traces of functions in W on these boundary parts. Indeed, 
w{x,t) = is in W, but its restriction to diK is not in K{diK). 

To study this further, dehne the maps 

Tj : n(x, t) H-)■ n(x, 0) and Tq : n(x, t) h-)■ n(x, x), 

whose application to any function gives its traces on d^K and doK, respectively. These 
maps are obviously well dehned for smooth functions. Below we prove that they extend 
to W. Let T^(S') denote the set of all measurable functions s on S' with hnite f^ws^. 

Lemma 4.1. For the W in this example, the following maps are continuous: 

Ti : IT-)■ L^(0,1), To : IT-)■ L^(0,1), and Ti - Tq : IT-)■ L^/^(0,1), 

i.e., there is a constant Cq > 0 such that 

/.I /.I /.I I |9 

\TiW — ToW\ 


(4.2) / 

2o 

for all w G IT. 


xItiwI dx + / xlrotcl ax + 


X 


dx < Collwll^ 


Proof. A general density result in [T51 Theorem 4] implies that C^{K) is dense in IT, so 
it suffices to prove (4.2) for all w G C^{K). Beginning with the fundamental theorem of 
calculus, 

nr 

Tpv{x) = w{x,r) — / dtw{x,s) ds, 

Jo 

squaring, integrating over r, and overestimating, 

nx nx nx nr 

x\tiw{x)\^ = / \tiw{x)\^ dr <2 / \w{x,r)\^ dr + 2 r \dtw{x, s)\‘^ ds dr. 

Jo Jo Jo Jo 

Now integrating over x and overestimating again. 


1 

2 


1 nx 


x|ritc(x)p dx < 


1 /•! 


\w{x, r)p dr dx + 


0 JO 


\dtw{x, s)|^ ds dr dx 


oJo 


= Uu 


w- 


A similar argument shows that the same inequality holds with r, replaced by Tq. 
To complete the proof, we therefore only need to show that 


(4.3) 


\TiW — ToW\ 


X 


dx < ||w||^. 
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But this follows from 

\TiW{x) — ToW{x)\‘^ = 


dtw{x, s) ds 


< X \dtw{x, s)p ds, 


dividing through by x and integrating over x. 


□ 


Lemma 4.2. Assumption (2.6) holds for this example after setting 

(4.4a) V = {w eW TiW = 0}, 

(4.4b) V* = {weW :ToW = 0}. 

Proof. For v,w ^ C^{K), the definition of D implies that 


1 px 


{Dw,v)w = 


{dtw)v + w{dtv) dt dx 


>0J0 

»i 


(4.5) 


{ToW){ToV) dx — {TiW){TiV) dx. 


In order to apply the density argument, we rewrite this expression: 


(4.6) {Dw,v)w= / (t^/^Tow) 


ToV — T\V 


X 


1/2 


dx + 


ToW — TiW 


X 


1/2 


{x^^'^TiV) dx. 


Now, one can immediately verify using Cauchy-Schwarz inequality and Lemma 4.1, that 
both the integrals extend continuously to W. Hence (4.6) holds for all v and w in W. 
Similarly, the expression 


(4.7) {Dw,v)w = 


>0 


ToW — TiW 


X 


1/2 


(x^^'^ToV) dx+ (x^^'^TiW) 


>0 


TqV — TiV 


X 


1/2 


dx, 


also holds for all v and w m.W. 


Let us verify (2.6a). For any v eV , since TiV = 0, we have from (4.7) that 


{Dv,v)w = 

Hence V C C~^. Similarly, V* C C~ 


ToV - 0 

xl/2 


(x^/^ToU) dx > 0. 


To prove (2.6b), let v eV. Then using (4.6) and putting TiV = 0, we have, 

/ /_'/; — ii\ 

dx 


{Dv*,v)w= I (x^^'^ToV*) ‘ 


>0 


X 


1/2 


which vanishes for any v* E V*. Hence V C '^D{V*). For the reverse inclusion, let 
G '^D{V*). Then, since ToV* = 0 for all v* eV*, we have from (4.6) that 

/■i 


{Dv*,v^)w = 


'0 


ToV — TiV 


X 


1/2 


(x^^^TiW^) dx = — iTiV*){TiV-^) dx. 

Jo 


Since all functions in D(0,1) can be written as TiV* for some v* E V*, this implies that 
TiW^ = 0 a.e. in (0,1), so E V. Thus, V = ~^D{V*). A similar argument shows that 
V* = ^D(V). □ 
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Remark 4.3. Note that although the two integrals in (4.5) need not generally exist for all 
w,v & W, those in the identities ( |4.6 ) and ( 4.7| ) exist for all G W. 

Remark 4.4. It is proved in [9l Lemma 4.4] that if 14 + 14* is closed, then an M that 

satisfies (2.4) can be constructed. They then write, “it is not yet clear to us whether 

properties ( 2.6a )-( 2.6b) actually imply that 14 + 14* is closed in W.” This issue was 
settled in [1] where they showed by a counterexample that (2.6a)-(2.6b) does not in 
general imply 14+14* is closed. Our study above provides another simpler counterexample: 
Specifically, for n > 2, let Xn denote the indicator function of the interval [1/n, 1]. Then 
Vn{x,t) = Xn{x)t/x is in 14 and n* = Xn{x){x — t)/x is in 14*. Clearly, as n —)■ oo, the 
sequence Vn + = Xn ^ V + V* converges in W. But its limit, the function 1, is not in 

14 + 14*. Indeed, if 1 were to equal v + v* for some v E V and n* G 14*, then by Lemma [4T| 


-dx = 

1 X 


WiV + ron*P 


dx < 2 


X 


\TiV\ 

X 


+ 


TnV 


* 12 


X 


dx < 2 C (||ii+ + |K'‘|Ih/ 


) 


which is impossible. 


Lemma 4.5. The inequalities of (2.11) hold for this example. 

Proof. Given any n G 14, by the density of C^{K) in IT, there is a sequence Wn G C^{K) 
converging to v in W. Let Vn{x,t) = Wn{x,t) — {T]Wn){x). Clearly G 14 fl C^{K). 
Moreover, 

\\xn — xWw < W^n — v\\w + ||||L2(K) 

( i-l px \ 1/2 

/ / \Wn{x, 0)P dtdx I 

= \\Wn — v\\w + ll+lL'nllLKO,!) 

= \\Wn — V\\w + \\Xi{Wn — 11 ) 1 ^ 2 ( 0 ^ 1 ). 


This, together with Lemma 4.1, and the convergence of Wn to v in W, imply the conver¬ 
gence of Vn to V in W. Thus 14 n C^{K) is dense in 14. Similarly, 14* fl C^{K) is dense in 
14*. Hence, it suffices to prove the inequalities of (2.11) for the dense subsets. 

For any C^{K) function v in 14, we have 


d 


d 


{x,tf = I —{v{x,sf) ds = I 2v{x,s)—v{x,s)ds, 


which implies 


1 px 


1 px 


v{x, tY dtdx <2 


oJo 


< 2 


>0 Jo 

f 

'0 


; \ 1/2 
v{x, sY ds I 


ds 


* d 

—n(x,s) 


X ft \ 1/2 

n(x, s)^ ds dt 


OJO 



\ 1/2 


ds ) dtdx, 




d , . 


1/2 


ds dt 


dx. 


Since x < 1, this shows that ||n|| 2 , < 2||Hn||i for all n G 14 fl C^{K) and hence for all 
V E V. The proof of (2.11b) is similar. □ 
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Theorem 4.6. Formulations (3.4) and (3.6) are well posed for this example 


Proof. By Lemmas 4.2 and 4.5, assumptions (2.6) and (2.11) hold, so Theorem |2.2| implies 
that (2.10) holds. Therefore, Theorem 3.1 gives the required result. □ 


4.2. Unidirectional advection. The above calculations have a straightforward gener¬ 
alization to multidimensional tent-shaped domains. We say that Ko is a vertex patch 
around a point p if it is an open polyhedron in (d > 1) that can be partitioned into a 
hnite number of d-simplices with a common vertex p G M'^. 

We hrst consider domains K built on (spatial) vertex patches of the form 


(4.8) 


K = [(x,t) : X e Ko, gi(x) < t < 9 „(i)} 


(and later, after Dehnition 4.8 below, specialize to tent-shaped domains). Above, go{x) 


and gfx) are Lipschitz functions on Kq such that iL is a nonempty open set in Then 

the unit outward normal vector n = exists a.e. on dK. Continuing to consider 


the same operator as in (4.1), namely A = dt, but on the new domain K, the following 


dehnes inflow, outflow, and characteristic parts of the boundary: 

(4.9a) diK = {{x,t) ^ dK : nt < 0}, doK = {(t, f) G dK : nt > 0}, 
(4.9b) dcK = {(x, t)edK :nt = 0}. 

We can immediately prove the following by extending the arguments of §4.1[ 


Theorem 4.7. Let K he as in (4.8) and let A = dt- Then the inflow and outflow trace 


maps, Ti : v{x,t) i—>• v{x,gi{x)) and Tq ■ v{x,t) i—)■ v{x, go{x)), extend to continuous linear 
operators 

rr-W ^ Ll_PKo), To-.W and r, - ^ 


As in (4.4), set V = ker(ri) and V* = ker(ro). Then, the assumptions of (2.6) and the 


inequalities of (2.11) hold. Hence the formulations (3.4) and (3.6) are well-posed. 


Identities similar to (4.6) and (4.7) prove the continuity properties of the trace maps 


stated above. To prove the stated wellposedness, we need to verify the assumptions 


in (2.6) and (2.11), which can be done by simple generalizations of the arguments in the 


proofs of Lemmas 4.2 and 4.5 Next, we proceed to consider a convection operator on 
tent-shaped domains. 


Definition 4.8. Suppose K and Kq are as in (4.8). If, in addition, K can be divided into 


hnitely many {d l)-simplices with a common edge {{p,t) : gfp) < t < go{p)}, then we 
call K a space-time tent. We refer to the common edge as its tent pole. Clearly, in this 
case, go and g^ are linear on each simplex of Kq. We split the tent’s boundary into the 
these parts: 


(4.10a) 

(4.10b) 


d;K = {(x, gfx)) : x G Kq], doK = {(x, go{x)) : x G Kq], 
d^K = dK\ {diK U doK). 
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We refer to the two parts in (4.10a) as the tent’s inflow and outflow boundaries, respec¬ 
tively. (Using such terms without regard to an underlying flow operator is an abuse of 
terminology that we overlook for expediency.) 

The equation modeling advection along a fixed direction a = (a,) G is of the form 
Au = f with 


(4.11) 


r\ d 

ou 




2 = 1 


du 

dxi 


Setting L = L‘^{K), D = T){K), and noting that A = —A, we can put this into the 


Friedrichs framework since the prerequisite (2.1) holds. 


Let n G denote the outward unit normal on dK. We often write it separating its 


space and time components as n = {n^^nt) with 
boundaries are such that 


We now assume that the tent 


(4.12a) diK C {(x, t) G OK : n at (x, t) satisfies n* -|- a • na, < 0} 

(4.12b) doK C {(x, t) G dK : n at (x, t) satisfies n* -|- a • na, > 0}. 

The vertical part of the boundary, namely c?biF, is further split into three parts d^K, 
and d^K where n^-l-acna; = a-n^ is positive, negative, and zero, respectively (see Figure]^. 
Let /] and /k denote the closures of d-^K U d^K and doK U d^K, respectively, and let 

r = r n r 

I lO - -ill I -i 0‘ 


Define 6{z) = dist( 2 ;, TJo). We will use the restriction of this function to /] and ik, as weight 
functions while describing the norm continuity of traces below. For smooth functions w 
on K, let 

nw = w\r„ ToW = w\r,. 


Theorem 4.9. Let K be a tent and A be given by (4.11). Suppose (4.12) holds. Then 
the above-defined maps t\ and Tq extend to continuous linear operators 

rr.W^L^fln) and To:W^L]{rfl. 


Hence V = ker(ri) and V* = ker(ro) are closed subspaces ofW. When restricted to these 
subspaces, the traces have an additional continuity property, namely 


(4.13) Ti : W ^ Ll/siT;) and r, : V ^ Ll/siTfl 


are continuous, 
well-posed. 


Finally, with this V and V*, 


the weak formulations (3.4) and (3.6) are 


Proof. The idea is to use a change of variable that brings the operator to the previously 
analyzed operator dt. The new variables are x = x — at and t = t, i.e.. 


X 

= H 

X 

t 


i 


where H 


I a 
0 1 
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Let K = H ^K. (Note that K is not a tent, in general.) Pnlling back fnnctions w on K 
to fnnctions w = w o H on K, the chain rnle gives 


(4.14) 


Aw = (Aw) o H, 


where A 


~dt 


Thns w e hh if and only A w & W = {z E L‘^{K) : Az G L‘^{K)}. 

Next, let n = {nx,nf) denote the unit outward normal on dK. Then n 
H^n/\\H^n\\ 2 . Dehning diK,doK, and dcK as in (4.9), we claim that 


(4.15a) dik = {{x,i) e dK : < 0} = H-\d;K U d^K) 

(4.15b) dok = {{x,i) e dk : fii > 0} = H-\doK U d+K), 

(4.15c) d^k = {{x,i) e dk : fii = 0} = H-Ad^K). 




■j 

O' 


V xgi 


V xgi 


1 _ 


_ -1 _ 


a ■ Vxgi - l_ 


For example, to sketch a proof of the hrst identity, note that n at {x,gi{x)) is in the 
direction of (V^S'i, —1) where Va, denotes the gradient with respect to x. Hence, because 
of (4.12), we have a ■ VxQi — 1 < 0 on d\k. Since the mapped normal n is in the direction 
of 

H^n = 

we conclude that < 0. Applying similar arguments on the remaining parts of the 
boundary, the claim (4.15) is proved. 

Let Kq be the projection of K on the t = 0 plane. There are (continuous piecewise 
linear) functions go and g; such that dok and dik a re gr aphs of go and gi, respectively, 
over Kq. On K, since A = df, we apply Theorem 
and foW = w 


to conclude that tiw = w 

extend to continous linear operators fj : IF —)■ L? (iFo) and f^ 


4.7 


diK 


W —)■ L?^_^.(iFo)- Hence V = ker(fi) and V* = ker('fo) are closed subspace of IF. By 


the additional continuity of A 
conclude that 


To : IF —)■ (also given by Theorem 4.7), we 


are also continuous. 

These continuity results are more conveniently mapped to K by using AA = dist( 2 ;, /jo). 
Note that go — gi vanishes at /jo = To restate the continuity properties of Tj in 

terms of 5, we prove that there are Ci, C 2 > 0 such that 

(4.16) ci(5(h,^i(h)) < ^o(h) -^i(h) < C2(5(h,^i(h)), Wx E Kq, 

(and similarly for Tq). When a point N = {x,gi{x)) on d^k is sufficiently near to /jo, the 
point P nearest to it on /jo, together with O = {x,go{x)) form a triangle (as shown in 
Figure]^. Now we may restrict ourselves to the two-dimensional plane containing this 
triangle. 

Consider the case when the segment PO lies on or below the plane of constant t passing 
through P, so that PO makes an angle 6(, > 0 with that plane. Let 6 be the angle made 
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K 


K 




Figure 2. On the left is a tent K with A = dt + 0.5dx that satishes (4.12) 


On the right is K obtained after applying the map in the proof of Theo- 
with mapped over operator A = df. 


rem 


4.9 


by PN and PO at P. Then, by elementary geometry, 

(4.17) 


7’-lV|l2=^l|0-iV|h, 

sint^ 


Note that 9 > 0 and 0 < 6*o < n/2. Therefore, observing that ||P — iV ||2 = S{N) and 
\\0 — A ^||2 = Qo — dh (4.17) proves ( |4.16 ). For the remaining geometrical conhgurations, 
identities similar to (4.17) can be derived to prove (4.16). Having established ( |4.16 ), we 
hnd that after mapping back to K, the stated continuity properties of TiW = {fiw) o 
and ToW = (fow) o H~^ are proved. 

It now onl y remains to prove the stated wellposedness of the weak formulations. By 
Theorem 4.7, for any F G hF', there is a unique u G L‘^{K) and q = Dz G (H*)-*- satisfying 


(4.18) - (h,ih)^2(^) + = F(i;), V h G IF, 

where D G T(IF, IF') is defined as before by {Dv,w)^ = {Av,w)j^2(^k) + L'^(k)- 


Here we have used (2.8) to hnd a i G F such that q = Dz. (While q is unique, i need 


not be unique.) It now follows from the properties of the mapping that u and q = Dz 
satishes (4.18[) if and only ii u = iio H~^ and z = z o satishes 


(4.19) 


{u,Av)l 2 (^k) + {Dz,v)w = F{v o H), V n G IF. 


Here we have used the fact that (4.14) implies (h, Hi))^ 2 (^) = {u, Av)l2(k) and conse¬ 
quently, {Dz,v)^ = {Dz,v)w This shows that the weak formulations on K and K are 


equivalent, so the wellposedness of the latter, namely (3.4), follows from the former. The 
wellposedness of (3.6) is proved similarly. □ 
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Remark 4.10. Under additional assumptions, including dist{diK,doK) > 0, a stronger 
trace result is proved in P Lemma 5.1]. However, on tents, dist(9iiL, doK) is always zero, 
so we are unable to use their result. 

4.3. A linear symmetric hyperbolic system in one space dimension. Let C G 

]^mxm Pe an m X m real symmetric matrix and let A be a tent as in 14.2 Set L = 
L^{K)^,T> = T>{K)^ and 


(4.20) 


, du ^ du 


where dtu and d^u are vectors in with their Uh component equal to dtUi and 


respectively. Since C is symmetric, A = —A, so assumption (2.1) is obviously satished. 


Let Q be an orthogonal matrix and A = diag(A£) be a diagonal matrix such that 
C = QAQ^. Let diK, doK and d^K be as defined in (4.10). In this subsection, we assume 


instead of (4.12) - that 


(4.21a) 

(4.21b) 


diK C {x G dK : ritl + UxC is negative definite}, 
doK C {x G OK : ritl + rixC is positive definite}. 


For each ^ = 1,... ,m, we decompose d^K into A, and d|]’^A where Xin^ is 

positive, negative, and zero, respectively. Then we have the following theorem, which is 
proved using the diagonalization of C to separate each component and then appealing to 


the analysis in § |4.2[ We now opt for a brief statement of the theorem, leaving the tacitly 
used properties of the traces to the proof. 


Theorem 4.11. Suppose (4.21) holds for the tent K and the operator A in (4.20). Then 


the formulations (3.4) and ( 3.6| ) with 
(4.22a) V = {z eW : = 0, for alii = 1,.. . ,m}, 

(4.22b) V* = {z eW : alli = l,..., m}, 

are well-posed. 

Proof. Let A = Q^AQ, W = {v G L‘^{K)'^ : Av G L^(A)”^}, and D be the corresponding 
boundary operator on W. Then clearly, n G IF if and only if h = Q^v is in W. Moreover, 
Aw = dtw + AdxW, i.e., its ith. component equals 

A(,wi = dtWi + Xedxiui. 


Note that Ai is a Friedrichs operator on A of the form (4.11) and has its associated graph 
space Wf and boundary operator D^. 


Now, the assumptions of (4.21) imply that (4.12) holds for each A^ (with a = Xi) 
so Theorem 4.9 yields the continuity of the maps r/ : tc h-)■ W(\pi and : w h-)■ Wi\ri 
on lU, where Tf = d^K U and Tf = doK U d^’^K. Therefore, the full trace maps 

Ti = (tj^, ..., r™) and Tq = (rQ,...,r™) are continuous on W. Set V = ker(ri) and 
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V* = ker(fo). Then the following variational equation for u G and q = Dz with 

^ e 

(4.23) - (m, Av)l 2 (k)-^ + {q, v)^ = F(h), Vh G W, 
splits into m decoupled equations, namely 

(4.24) - {ui, At>Vi)L2,^K) + {DeZi, Vi)^^ = F{ve), Wve eWe, = 



To transfer these results for A to A, we dehne 


TiW = fi{Q^w), ToW = foiQ^w). 


Then (4.22) is the same as Id = ker(ri) and V* = ker(ro). Note that ^ G Id if and only if 
z = Q^z G Id. Also note that a u G and ^ G Id solves (4.23) if and only ii u = Qu 

and 2 ; = Qz satishes 

-(m, An)i 2 (;^)m + {Dz,v)w = F{Q^v), Vn G W. 

Here we have used {u, Av)l 2 (^j^^ m = {u,Av)l2(k) m and consequent identities for the cor¬ 
responding boundary operators. Thus the stated wellposedness of ( |3.4 ) follows from the 
established wellposedness of (4.23). The proof of wellposedness of (3.6) is similar. □ 


Remark 4.12. Consider a tent K with empty d\^K. Then, under the assumptions of 
Theorem 4.11 , a function in Id has all its m components equal to zero on the inflow 
boundary diK. Moreover, if n G Id fl C{K), then applying the additional continuity 


property (4.13) to the operators in the above proof, we hnd that the outflow trace of 


each component of v must approach zero as we approach Ro where the inflow and outflow 
boundary parts meet. 


5. The wave equation 


We now apply the previous ideas to the important example of the wave equation and 
work out the resulting weak formulation in detail. Our model problem is to hnd a real¬ 
valued function 0 on the space-time domain i7 = (0,5') x (0,T), satisfying 


(5.1a) 

c dttcf) dxx4^ S') 

(5.1b) 

o' 

(5.1c) 

dt(j) - cdcccj) = 0, 

(5.Id) 

dt(j) + cdx(j) = 0, 


0<a;<5, 0<f<T, 
f = 0, 0 < a; < S', 
a; = 0, 0 < f < T, 

X = S, 0 < t < T, 


where c > 0 is the wave speed. Here, we have imposed the outgoing impedance boundary 
conditions (but other boundary conditions can also be considered - see Section]^. 
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The above second order system for 0 arises from a system of first order physical princi¬ 
ples, which also matches the form of the problems we have been studying, namely (2.9). 
Set 

cd^cj) 


u = 


dtcf) 


and observe that dtUi = cdxt4> = cdxU 2 and dtU 2 = dtt<p = cdxUi + c^g. These two 
equations give the first order system Au = f where 


(5.2) 


Au = dtu 


0 c 
c 0 


dxU, f = 


0 

c^g 


It hts into the framework of 14.3 after the diagonalization 


0 c 
c 0 


= QAQ\ Q = ^ 


V2 [1 


A = 


Ai 0 
0 Aa 


C = - 

where Ai = —c and Aa = c. 

Analogous to (4.9), we define diQ = (0, S') x {0}, dof2 = (0,5') x {T} and = 
df2 \ {dif2 U c?of?)- The vertical parts are further split into 


= {0} X [0,T], = {S} X [0,T]. 






5+i2 I 


Set rf and to the closures of and respectively, Ff^ = FfdF^, 

and = dist((a;, t), i^f,) for ^ = 1,2. By a minor modihcation of the arguments in 

Section one can prove that the global trace maps 


Zi 


[Zi + Z2)\pi 

A2. 


1 

1 

0 

to 

. 

to 

1_ 



{zi + X2)| 

^ 0 

A2. 


1 

1 

0 

to 

Oto 

1_ 


are continuous. Set 


rr.W^ Ll{Fl) X Ll{F^) 
To : IT ^ Ll{Fl) X Ll{Fl) 


V{Q) = ker(ri), V*{Q) = ker(ro). 


These spaces can be used to give a global weak formulation on i7, but our focus in on 
local solvers. 

In space-time tent pitching methods, we are required to numerically solve the wave 
equation on space-time tents, ordered so that inflow data on a tent can be provided by 
the outflow solution on previously handled tents or through given data. Hence we now 
focus on the formulation and discretization on one tent K. 


5.1. Weak formulation on a tent. Consider the analogue of (5.1) on one tent K, with 
zero initial data on the inflow boundaries and with boundary conditions inherited from 
the global boundary conditions (5.1c)-(5.1d). 
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Define, as before, the boundary parts of a tent K, by 
diK = {(x, t)edK :nt< 0}, doK = {(x, t) e OK : ut > 0}, d^K = dK\ {d^K U doK), 

d^’^K = d-'^K = {(x,f) e d^K : cn, < 0}, 

d^’^K = = {{x,t) G dijK : cn^ > 0}. 

Note that the boundary part d\^K may be empty in some tents. We consider the tent 
problem of solving for u satisfying 


Au = f 
n = 0 


on K, 
on 


Ml — ^2 = 0 
Ml + M2 = 0 


on 

on 


To obtain a well-posed weak formulation on one tent, we proceed to use Theorem 4.11 


To this end, we must assume that the tent satisfies Assumption (4.21), which now reads 


(5.3a) 

(5.3b) 


d-^K C {x G dK : rit ± n^c < 0}, 
doK C {x G dK : rit ± n^c > 0}. 


Since A = —A in this example, the weak formulation (3.4) reads 

(5.4) u e L, q e : - {u, Av)l + {q,v)w = K{v), Vm G IT, 


where the spaces are set following (4.22), namely 



^1 

/2. 

G W : 

Zl + Zi 

_zi - Z‘i_ 

£ 

T* = 1 

Zl 

A2. 

G W : 

Zl + Z2 
_Zi - Z2_ 

£ 


= 0, for £=1,2 


diKud~’‘K 


= 0, for £=1,2 


(9oAuak A 


Theorem 4.11 shows that (5.4) is a well-posed weak formulation on K provided the tent 


K satishes (5.3). Note that the above spaces change from tent to tent and may arguably 


be better denoted by V{K),V*{K), etc., but to avoid notational bulk we will suppress 
the A'-dependence. 


5.2. CFL condition. Let us take a closer look at (5.3). First note that each tent, in this 


application, consists of either two triangles (on either side of the tent pole), or just one 
triangle. The tents are thus divided into three types, as shown in Figure 

The length of the tent pole is /c, the numbers pi and pr are such that pik and prk give the 
heights of the outflow boundaries on the left and right side of the tent pole, respectively, 
and the spatial mesh size are hr, hi > 0. Writing down the normal vector on the tent 


boundaries, we immediately hnd that condition (5.3) on a tent is equivalent to 
(5.5) c -— < 1 and c—— < 1. 

hr hi 

Clearly, by controlling the size of the tent pole we can satisfy these inequalities. 
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Type I: > 0, /i/ > 0 


Type L: > 0, /i; = 0 Type R: hr = 


0, /i/ > 0 


Figure 3. Three types of tents 


The well-known Courant-Friedrichs-Levy (CFL) condition [5] identifies stability con¬ 
ditions as constraints on the time step size in terms of space mesh size in nnmerical 
discretizations. In onr case, this condition manifests itself as geometrical constraints (5.5) 


on the tent. For this reason, we will refer to (5.3) 
onr method. 


or (5.5) - as the CFL condition of 


6. The numerical scheme 


In this section, continning to work with the wave operator A defined by (5.2), we give 
an explicit nnmerical scheme for approximating u{x,t) satisfying 


(6.1a) Au = /, 

(6.1b) Ui{x,0) = u'^{x), 

(6.1c) Mi(0,t) - M 2 ( 0 ,t) = 0, 

(6.Id) Mi(l,f)-f M2(1,^) = 0, 


0<x<5, 0<f<T, 
0 < a: < S', i = 1,2, 
0<t<T, 

0<t<T. 


The scheme will allow varying spatial and temporal mesh sizes. Here / and are assnmed 
to given smooth fnnctions. We begin by describing the calcnlations within each tent, 
followed by the tent pitching techniqne to advance in time. 


6.1. Conforming discretization on a tent. As seen above, a tent is comprised of one 
or two triangles. Let the space of continnons fnnctions on a tent K whose restrictions to 
these triangles are linear be denoted by P^{K). We construct a conforming discretization 
of ( 5.4[ ) within K using the discrete space 

(6.2) Hi = F n 


By definition, Vi C V, and consequently, functions in Vi must satisfy the essential bound¬ 
ary conditions of V. Depending on the tent geometry, different boundary conditions must 
be imposed on different tents. 

To examine what this entails for the nodal coefficients on mesh vertices, let ( G Pi{K) 
be the continuous scalar function (unique Lagrange basis function) that equals one at the 
“apex” of the tent K, equals zero at all its other vertices. The apex of a tent, irrespective 





















22 


GOPALAKRISHNAN, MONK, AND SEPULVEDA 


of whether it consists of one or two triangles, is the vertex in doK that is away from diK. 


Now, snppose fi = 


— [Ml ■ 
.M2 . 


in 


is snch that 


if d^’^K is nonempty. 


(6.3a) /xi - /r 2 = 0 

(6.3b) /ii + = 0 if d^’^K is nonempty. 

(Note that if d^K is empty, then /r is an arbitrary vector in M^.) Then, it is easy to see 
that 

(6.4) Vi = : /i satishes (6.3)} 

provides an alternate characterization of ( |6.2 ). 

A compntable conforming discretization of (5.4) additionally requires hnite-dimensional 
subspaces of L and (V^*) . For the latter, observe that (2.8) implies that 

D{Vi) C D{V) = 


Hence we choose an approximation qi of the solution component q in (5.4) to have the 
form 

qi = Dzi, Zi E Hi. 

Then qi is clearly in (H*)"''. Next, set Li C T to be the space of vector functions whose 
components are constant functions on K. Finally, set 

(6.5) Wi = [w-. tc = K + pC, A, k(^EVi]. 

Our discretization of (5.4) now takes the following form: Find Ui in Li and qi E T)(Vi) 
satisfying —{ui,Aw) + {qi,w) = F{w), for all w E Wi. Clearly, dim(hFi) is four or three, 
depending on whether d\,K is empty or not. This equation gives rise to an invertible 
discrete system, as a consequence of the unisolvency of the following slightly modified 
problem: 

Find ui E Li and Zi E Vi such that 

( 6 . 6 ) 

— {ui, Aw) + {Dzi, w) = F{w), Vtc G Wi- 


Proposition 6.1. There is a unique solution for Problem (6.6). 


Proof. Note the dim(Li) +dim(Vi) = dim(hFi), so (6.6) gives a square (Petrov-Galerkin) 
system. Hence it suffices to set F = 0 and prove that Ui = Zi = 0. With F = 0, writing 
Zi = af for some a G we have {D{a(), w) = 0 for all constant w E W. Since aC, eVi, 
by the definition of IFi, we may set w = a in (6.6), to get 

{D{aQ,a) = [ ■ aC = 0 

where 

(6.7) ^ = Ufl + Cux = 


'doKudi^K 


nt —cux 
—CHx nt 


If d\,K is empty, then since ■ a = (oi + a 2 )^{nt — cnx)/2 + (oi — Q; 2 )^(^t + cnx)/2, the 
CFL condition (5.3) gives a = 0. If d^K is nonempty, then whenever G H we have 
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either cii — 0^2 = 0 or cii + 0^2 = 0, so we can continue to conclude that a = 0. Of course 
a = 0 implies Zi = 0. 

To prove that Ui = 0, we use ( 2.3| ) after substituting zi = 0, to get 

{Dw,ui) = 0, \/w G Wi- 

Since ui is a constant function, uiC G Wi, so we may choose w = uiC and conclude that 
Ml = 0 by an argument analogous to what we used above. □ 

Remark 6.2. One can view zi\d^K as an interface trace variable and qi = Dzi as an 
interface flux variable. By the trace theory we developed previously, outflow trace zi\q^k 
must vanish at the points where outflow and inflow edges meet in order for zi to be in 
V. This motivates our choice (6.4) of Vi to obtain a conforming method. Other non¬ 


conforming avenues to design approximations within a tent can be found in [10] and [19] . 

6.2. Advancing in time by tent pitching. We now show how the above ideas yield an 
explicit time marching algorithm for solving (6.1). First, we mesh the space-time domain 
i7 = (0, S) X (0,T) by a collection of tents K with these properties: The hrst property 
is that either d\^K is empty or 

(6.8a) dbiF C 5bi?, 

for all K G Qh- Second, there exists an enumeration of all tents, iFi, A 2 ,... ,iFj, with 
the property that for each j G {1,... J}, 

(6.8b) d^Kj C IJ doKk U d^n. 


k=l 


Finally, for all j G {1,... J}, 


(6.8c) 


Kj satishes the CFL condition (5.5). 


It is well-known how to construct an algorithm (not only in one space dimension, but also 
in higher dimensions [SI 1^ ) that produces meshes satisfying ( |6.8D , so we shall not dwell 
further on the meshing process. 

The discrete space-time approximation on the mesh is developed using 
(6.9a) 14 = {x G n V{f2) : z\k G VA G A/,, and z{x,0) = hu^x)} 

(6.9b) Lh = {a : a\K G is constant on each K G f2h} 

(6.9c) Wh = {w : w\k & IFi on each K G i7/i} , 

where denote the linear nodal interpolant on the spatial mesh. The method hnds 
approximations Uh G Lh and Zh E Vh satisfying 


( 6 . 10 ) 


^ ( - / Uh- Aw + &Zh ■ / 


f -W, VtC G Wh, 


K 


where ^ is as in (6.7). 

Because of (6.8), we are able to use a time-marching algorithm to solve (6.10): Proceed 
in the ordering of (6.8b), and for each tent K, solve for Uh\K and Zh\K- Specihcally, if a 
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is the nodal (vector) value of Zh at the apex of K, then dehning the problem on 

one tent is to hnd Uhlx £ -hi and z^ G Vi satisfying (6.10), namely 


( 6 . 11 ) 


Uh ■ Aw + 


IK 


/ ^z^ ■w= / f -w- 

'dK Jk JdK 


^zi 


K 


W, 


yw e Wi. 


where zf^ = Zh — z^. Note that zA on right hand side will be a known quantity if (6.8b) 


holds and if we have already solved on every K' appearing before K in the ordering of 
tents in (6.8). Indeed, z^ is completely determined by its nodal values at (the three or 


two) vertices on diK, which either lie at f = 0 or were apex vertices of previous tents. 


Problem (6.11) is exactly of the same type we discussed in § 6.1 


6.3. Propagation formula. Since the system (6.11) is small, we can explicitly calculate 


its solution. To see how information is propagated from inflow to outflow on a mesh of 
tents, we consider the case where the volume source / is zero. Write Zh = in (6.10) 

and let the nodal values of the scalar Lagrange hnite element functions Zh^i and Zh ,2 be 
[yl~\ , [yl] , [yi] , [yr], at the top, bottom, left and right vertices, respectively, of a tent 
of Type I, as in Figure For the other two tent types, we omit the nodal values at the 
missing vertex. 

Equation ( |6.11| ) hnds [yl] as a function of the remaining nodal values. After tedious 
simplihcations (not displayed), this relationship is found to be as follows: 


for Type I, 
for Type L, 
for Type R, 

for Type I, 
for Type L, 
for Type R. 







0 c' 




'U^ - U^' 

(6.12) 

_w_ 

— 

y\ 

+ Wi 

c 0 


- v‘ 

+ W2C 

- v‘ 






0 c' 


■f/r- _ l/b- 


■f/r- _ l/b- 


y\ 

— 

y\ 

+ Wi 

c 0 


yr _ yb 

+ W2C 

yr _ yb 






0 c' 


'U^ - [/'■ 


- y 


_w_ 

— 

y\ 

+ Wi 

c 0 



+ W2C 

yb - y^ 


where 

Wi = 

Wi = 

Wi = 


{hr + hi)k 


{hr + hiY - c^k‘^{pr - piY 

k 

2{ck{l — Pr) + hr) ’ 

k 

2{ck{l -pi) + hi)' 


, W2 = 


Ck‘^{pr - Pi) 


{hr + hz)2 - C^k‘^{pr - Pl)"^' 
W2 = Wi, 

W2 = -Wi, 


6.4. Error analysis on uniform grids. We now work out the stencil given by the 
method on a uniform grid where all tents are shaped the same (see Figure]^. The stencil 
translates (6.12) into an equation that gives the nodal values of the outflow apex vertex. 


given the nodal values at the inflow vertices. Let h > 0 be the uniform spatial mesh size, 
fc > 0 be the time step size measured, as before, by the height of the tent pole. At a 
point {jh/2, kn/2) in the lattice {h/2)'L x {k/2)'L, let (f/”. Id") denote the nodal value of 
the approximation to Zh there. As shown in Figure the scheme uses only a subset of 
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A ;/2 


Figure 4. The stencil 


lattice points in KL x fcZ. Each grid point involved in the scheme has an associated U 
value (indicated in the hgure by “•”) and a V value (indicated by “o”). 


Equation (6.12) now simplihes to 
(6.13a) 

(6.13b) 


^n+l ^ jjn-l ^ 


v; 


n+1 


= 1/’ 


n—1 


/+! - U-.) 

+ ac(f/yi - 


where a = k/h. This is simply the non-staggered leapfrog scheme (studied extensively for 
scalar equations) applied to the hrst order system. By a simple Taylor expansion about 
the stencil center, we see that the scheme is consistent and that the local truncation error 
is of second order (see |23l [21] for dehnitions of these and related terminology). 


To examine stability, introduce a new vector variable X"’ and rewrite the scheme (6.13) 


as follows: 




\un 

|A'”]4 + ocl.y+i - 

, where X" = 

yn 

[A7]i 


Un-1 

[ATI. J 


Wri 


. The amplihcation 


xr+i = 


To this one-step scheme, we now apply von Neumann analysis 
matrix G, connecting to X” can be readily calculated: 


G = 


% denotes the imaginary unit, the eigenvalues of G are g\ = is — a/1 — g 2 = -\- 


■ 0 

2is 

1 

O' 



■ 1 

1 

1 

1 ■ 

2is 

0 

0 

1 

= RAR-\ 

where R = 

1 

1 

1 

1 

1 

0 

0 

0 





. 0 

1 

0 

0. 



.^1 ^ 

92^ 


-94 \ 


yr 


Sb 5'3 = - vT- s% gi = -IS 

9 G [—vr, tt] gives the frequency in von Neumann analysis. If 


-f a/1 — ty = diag( 5 fj), s = acsin(6*), and 


(6.14) 


lacl < 1 
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h-convergence 




Number of spatial mesh points (1/h) 
(b) Second order convergence 


(a) Errors as a function of time 

Figure 5. Comparison with CTCS scheme 




then all eigenvalues satisfy \gi\ = 1. Furthermore, since det R = A{g^ — g 2 )((y '4 
remains away from zero whenever (6.14) holds, the powers G" = RA^R~^ are uniformly 


bounded for all n and all 9. Hence (6.14) implies that the scheme is stable. 


We thus conclude, by the Lax-Richtmyer theorem, that the scheme is convergent and 
is of second order. Note that the CFL condition we previously found on general meshes. 


namely (5.5), when restricted to uniform meshes, gives exactly the same CFL condi¬ 


tion (6.14) obtained above from von Neumann analysis. 


7. Numerical results 


7.1. Convergence study. First, we report numerical results from our tent pitching (TP) 
scheme and compare it with the well-known “central-time central-space” (CTCS) finite 
difference scheme (see [5l [23], sometimes also known as the Yee scheme [121 EZ])- The 
only difference between the two is that while the TP scheme sets the U and V nodes 
on the same location (exactly as indicated in Figure]^, the CTCS scheme sets them 


on 


staggered locations on the same grid. Both schemes are applied to the model problem (6.1) 
on uniform grids with S' = 1. We use a grid like that in Figure]^ for both methods. 

To impose the outgoing impedance boundary conditions within the CTCS scheme, we 
use the standard finite difference technique of introducing ghost points to the left and 
right of the hnite grid and eliminating the unknown values at those points using the 
boundary condition. In contrast, in the TP scheme, the impedance boundary conditions 
are essentially imposed within the hnite element spaces, as we have already seen previously. 
We set c = 1 and impose the initial condition so that the exact solution is 

ui(a;,t) = U2(a;,t) = 


i.e., the solution is a smooth pulse moving to the left at unit speed, eventually clearing 
out of the simulation domain. At every other time step (in the uniform space-time grid) 
we compute the T^(0, l)-norm of the difference between the computed and exact solution. 





























SPACE-TIME TENT PITCHING AND TRACES 


27 


The evolution of these errors in time on a grid of spatial mesh size h = 0.0025 and k = 0.9h 
is shown in Figure |5^ 

We observe from Figure [5^ that the errors of both methods are comparable and remain 
low throughout the simulated time. Note also that after the pulse clears the simulation 
domain reflectionlessly (and the solution within [0,1] vanishes), the errors for both meth¬ 
ods decrease markedly. In Figure we display a log-log plot of the T^(0, l)-norm of the 
errors at f = 0.5 for h = 1/2^,..., 1/2^^, and k = 0.9h. The rate of decrease of this error 
is clearly seen to be of the order O(h^). This is in accordance with our von Neumann 
analysis of § 6.4 (although we did not take into account boundary conditions in that 


analysis). 

Thus we conclude from Figure that there is negligible difference between the perfor¬ 
mance of the two methods on uniform grids. 


7.2. Material interfaces and other bonndary conditions. Next, we consider a gen¬ 


eralization of (6.1) given by 


(7.1a) 

(7.1b) 

(7.1c) 

(7.1d) 

(7.1e) 


dt 

Ki 

o' 


Ui 


' O C 

o 

dx 

Ml 


0 

^2 


U2 


c U 


U2 


ui{x,0) = ul{x), 
U2{x, 0 ) = ul{x), 
ZqUi — U2 = 0, 

ZiUi +U2 = 0, 


0 < X < 1,0 < t < T, 
0 < X < 1, 

0 < X < 1, 

X = 0, 0 < f < T, 

X = 1, 0 < f < T. 


where ki{x) and K 2 {x) are time-independent material parameters and c, zq and zi are 
constants. Such systems arise from electromagnetics or acoustics [121 US] on layered 
media and the differential equation is often written in the following equivalent, but non- 
symmetric form 


Ml 


'0 

A' 

dx 

Ml 

“2. 


P 2 

0 


“2. 


/ 


where (3i{x) = c/Ki{x) and / = diag(K7^, f^ 2 ^)f obtained by scaling the equations of (7.1a) 
by ^7^ and When Ai(x) = H 2 {x) = 1 and Zq = Zi = 1, we obtain the model formu¬ 
lation we discussed previously in detail. Dirichlet boundary conditions can be imposed 
by putting zq = Zi = 0, while exact outgoing impedance conditions can be imposed us¬ 
ing zq = ^ki/k 2 and Zi = \/ni/ k, 2 - Intermediate values of Zi give damped impedance 
boundary conditions. 

Whenever Ki is a constant on each spatial mesh interval, a tent pitching scheme is 
suggested by a simple generalization of the previous algorithm for homogeneous media. 
We define the discrete spaces exactly as in (6.9), but noting that V{Q) now has different 
essential boundary conditions - stemming from (7.1d)~-(7.1e) - which are inherited by the 
spaces on tents with its tent pole on the boundary. The generalization of the scheme is 
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at timeO.OOO 



X 


(a) Time snapshots of ui (cf. whole 
space-time plot of ui in Figure [6b|) 



(b) Entire space-time plot (over all time slabs) 
of the ui-component of the solution 



(c) Parts of the tent pitched mesh of the initial 
time slab, with varying spatial mesh sizes in the 
regions x < 0.5 and x > 0.5. 


Figure 6. Wave propagation through an impedance-matched interface 


derived by merely setting the A and ^ in (6.11) by 


A = 


Ki 

0 


0 

^2 


a 


a 


ritKi —CHx 
— CTlx ntK2 


Note that this A, appearing on the left hand side of (7.1a), satishes (2.1). By solving this 
general version of (6.11) one can obtain propagation formulas similar to (6.12), but we 
omit these details and report only the numerical results. 

First we consider the case 


Kl = 


K2 = 


0 < X < 1/2, 
1/2 < X < 1, 


2, 0 < X < 1/2, 

1, 1/2 < X < 1, 

and c = 1. The wave speed (equalling c/yTTiF^), jumps from 0.5 in the left half to 1 in 
the right half. However, the impedance (equalling ki/K 2 - see [I6]) is one in both regions. 
Thus X = 0.5 is an impedance-matched interface about which we do not expect to see any 
reflection. 

We use the tent pitching method to simulate a wave propagating to the right starting 
near x = 0.2. To this end, dehne a smooth pulse g{x) = and set the data 

in (lO by 


(7.2) 


/ = 0, M?(x) = {c/Ki)g{x), ul{x) = -{c/y/K^)g{x), 


and zq = y/ ki/k 2 and zi = y/ ki/k 2 . We use a spatial mesh of mesh size h = 10 ^ in 
the left half and h = 2 x 10“^ in the right half. A simple tent meshing algorithm then 
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produces a mesh of space-time tents based on this non-uniform spatial mesh that satisfies 
the CFL condition (5.5). The meshing algorithm proceeds as illustrated as in Figure 
by simply picking a point with the lowest time coordinate to pitch a tent. When multiple 
locations have the minimal time coordinate, the algorithm picks a tent pitching location 
among them randomly, thus giving an unstructured mesh. To minimize the overhead in 
constructing the mesh of tents, instead of meshing the entire space-time domain at once, 
we hrst mesh a thin time slab {(x, t) : 0 < t < 0.002, 0 < x < 1} and then repeatedly 
stack this mesh in time to cover the entire region of time simulation. The mesh of the 


initial slab is shown in Figure 6c 


One of the two components of the computed solution is shown in the remaining two 
plots of Figure Clearly, the simulated wave packet travels left across the x = 0.5 
interface without any reflected wave and expands as it enters the region of higher wave 
speed. In further (unreported) numerical experiments, we have noticed changes in the 
discrete wave speed depending on the space-time mesh. For example, the wave speed 
differs if one uses uniform space time meshes with positively sloped diagonals only or 
negatively sloped diagonals only. Such wave speed differences appear to approach to zero 
slowly as h is made smaller. High order methods may be needed to reduce these dispersive 
errors. 

Our next and final example involves an interface where we expect both reflection and 
transmission. We set c = 1 and 


Ki = 


K2 = 


1 , 

1 / 2 , 


0 < X < 1/2, 

1/2 < X < 1, 


4, 0 < X < 1/2, 

1/2, l/2<x<l. 

Both the wave speed and the impedance jumps from the left region to the right region 
(from 0.5 and 4 to 2 and 1, respectively). We set / and initial data as in the last simulation 
by dz but in order to impose Dirichlet boundary condition, we set zq = zi = 0. This 
time, instead of using a non-uniform mesh, we use a spatially uniform mesh of h = 10“^ 
and let the tent pitching algorithm adjust k to satisfy the CFL condition (5.5) in each 
tent. We found that the mesh obtained, displayed in Figure 7a, while not ideal due to 
the thin triangles, is adequate for the simulation. (Better tent pitched meshes can be 
obtained using non-uniform spatial mesh spacing, as we saw in the previous example 
The solution components Ui and U 2 obtained from the simulation are 
The computed waves are transmitted as well as reflected both 


and Figure 6c 


displayed in Figure 7b 


from the interface and the Dirichlet boundaries. The expected features of the solution are 
therefore recovered by the method. 
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